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Chapter  1 


Introduction 


With  the  prospects  of  space  travel  so  near,  evaluation  and  design  of  structures  suit¬ 
able  for  a  space  environment  is  needed.  These  structures  include  satellites,  the  space 
station,  space  antennas,  and  others,  and  are  grouped  in  a  category  termed  Large 
Flexible  Space  Structures  (LFSS).  Since  many  LFSS  are  characterized  by  low  damp¬ 
ing,  vibration  control  is  required  to  achieve  sufficiently  accurate  shape  control.  There 
are  many  ways  to  approach  vibration  control,  including  passive  damping  and  active 
damping.  In  the  case  of  active  damping  through  modal  control,  specific  modes  which 
contribute  significantly  to  the  unwanted  motion  are  selectively  damped  through  the 
use  of  active  controls.  A  sensor  senses  the  motion  and  processes  it  through  a  filter, 
and  a  controller  sends  the  signal  to  an  actuator  in  an  effort  to  counterbalance  the 
undesired  modes. 

The  difficulties  encountered  in  the  evaluation  and  design  of  LFSS  control  systems 
stem  from  the  attempts  to  accurately  model  such  a  structure.  Using  modal  control 
involves  formulating  the  model  in  modal  form.  Technically,  structural  modes  are 
a  mathematical  representation  of  actual  motion.  The  motion  is  actually  composed 
of  an  infinite  number  of  modes,  for  which  controllers  are  impossible  to  design  or 
to  implement.  Furthermore,  retaining  a  large  number  of  modes  is  not  practical  for 
computational  purposes.  Nevertheless,  a  “truth”  model  of  the  LFSS  must  contain 
a  reasonable  number  of  modes  for  accurate  modelling  of  the  dynamics  and  inherent 
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characteristics  of  the  LFSS.  This  modelling  process  is  easier  said  than  done  -  these 
structures  are  difficult  to  test  on  earth  due  to  their  size  and  the  lack  of  a  suitable 
space  environment.  The  structural  frequencies  of  LFSS  are  generally  low,  hence 
easily  excited,  and  closely  spaced  which  creates  difficulties  in  modelling  due  to  the 
inability  to  distinguish  the  modes  of  the  structure.  Many  experiments  [1-3]  are  being 
executed  on  scaled-down  models,  some  with  zero  gravity  springs,  to  determine  the 
characteristics  of  these  structures  -  mainly,  frequencies  of  the  structural  modes,  mode 
shapes,  and  damping  of  the  structure.  Since  these  experiments  contain  errors,  as  well 
as  the  fact  that  the  dynamics  and  characteristics  of  the  structure  may  change  with 
time,  there  are  uncertainties  in  the  parameters  of  these  structures.  Therefore,  the 
truth  model  of  the  LFSS  may  not  adequately  represent  the  structure  through  all 
operational  conditions. 

One  way  to  improve  the  truth  model  is  to  introduce  adaptive  estimation,  which 
modifies  the  truth  model  through  actual  information  obtained.  Unfortunately,  con¬ 
vergence  to  a  solution  is  not  always  guaranteed.  Also,  adaptive  estimation  tends  to 
be  more  complex  and  increases  computational  difficulties. 

Another  method  is  to  design  a  control  system  to  be  “robust”  to  modelling  uncer¬ 
tainties.  In  this  manner,  the  closed-loop  system  will  be  asymptotically  stable  in  the 
presence  of  these  uncertainties.  Two  major  problems  axe  encountered  with  robust 
designs.  The  first  is  that  so  many  assumptions  are  made  that  the  reality  of  the  LFSS 
hardware  is  evaded  and  the  control  system,  under  these  circumstances,  may  not  be 
robust.  It  is  assumed  in  this  study  that  the  actuators  and  sensors  are  point  masses, 
velocity  sensors  are  employed,  and  the  dynamics  of  the  actuators  and  sensors  axe 
much  faster  than  those  of  the  dominant  modes  of  the  structure  itself,  allowing  the 
higher  frequencies  to  be  filtered  out.  Also,  all  actuators  and  sensors  are  collocated 
-  in  this  way,  the  force  can  be  applied  directly  to  the  point  where  the  measurement 
of  motion  was  taken.  Mathematically,  this  collocation  ensures  the  plant  “positive 
realness”  and,  hence,  desirable  stability  properties  [4].  One  of  these  properties  states 
that  for  a  positive  real  plant,  if  the  feedback  design  is  strictly  positive  real,  then  the 
closed-loop  system  is  robust.  Thus,  although  the  actuator  and  sensor  pairs  or  groups 
of  pairs  may  have  mass  and  may  have  “lower”  frequency  dynamics,  it  is  hoped  that 
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the  effects  get  washed  out  with  the  guarantees  of  robustness.  Furthermore,  the  ad¬ 
vancement  of  technology  in  space  actuators  and  sensors  will  hopefully  reduce  many 
concerns.  The  second  limitation  with  robust  designs  is  that  the  control  system  may  be 
designed  to  be  so  insensitive  to  parameter  changes  that  performance  is  not  optimal. 
There  is  indeed  a  tradeoff  between  sensitivity  and  performance,  which  is  addressed 
in  Chapter  4. 

To  implement  modal  control  [5],  estimates  of  the  modal  positions  and  velocities 
of  the  structure  are  required.  Full  and  reduced-order  centralized  Kalman  filters  have 
been  reviewed  extensively  in  [6-Sj.  The  stability  and  robustness  of  decentralized 
filters  is  discussed  in  [9-  14].  These  studies  do  not  emphasize  Single-Input  Multi- 
Output  (SIMO)  systems  which  are  a  subset  of  Multi-Input  Multi-Output  (MIMO) 
systems.  This  study  focuses  on  the  performance  of  optimal  and  suboptimal  filters 
under  such  circumstances  with  perfect  and  imperfect  structural  frequency  informa¬ 
tion.  The  estimators  evaluated  here  include  the  optimal  Kalman  filter,  which  is  a 
full-order  centralized  filter,  centralized  reduced-  order  Kalman  filters,  a  decentralized 
filter  in  which  error  amplitude  frequency-shaping  occurs,  and  sensitivity-shaped  fil¬ 
ters  designed  to  reduce  sensitivities  to  variations  in  the  structural  frequencies.  This 
investigation  includes  systems  with  and  without  modelled  noise. 

Chapter  2  introduces  optimal  and  suboptimal  estimators.  Chapter  3  contains 
modal  analysis,  problem  formulation,  and  perfect  model  information  results  for  cen¬ 
tralized  and  decentralized  estimators.  Section  4.1  proposes  one-mode  and  two-mode 
sensitivity-shaped  filters  which  are  less  sensitive  to  errors  in  the  structural  frequen¬ 
cies.  The  sensitivity-shaped  filter  performance  is  compared  with  the  centralized  and 
decentralized  estimators  of  Chapter  3  under  structured  frequency  variations  of  20  %. 
Section  4.2  presents  basic  concepts  of  positivity  [4].  The  robustness  of  the  system 
under  structural  frequency  errors  with  the  one-mode  sensitivity-shaped  estimator  is 
analyzed  using  a  phase-shaping  method  of  Section  4.3  and  the  Lyapunov  method  of 
Section  4.4.  The  phase-shaping  method  is  based  on  the  principle  of  positivity  [4]  with 
collocated  actuators  and  sensors.  It  is  a  conservative  method  which  requires  the  phase 
of  the  feedback  loop  estimator  transfer  functions  to  lie  between  ±  90  degrees,  along 
with  conditions  on  the  controller  gains,  for  the  closed-loop  system  to  be  asvmptoti- 
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cally  stable.  The  sensitivity-shaped  filter  is  shown  to  be  at  least  marginally  stable  by 
this  method.  The  Lyapunov  method  requires  the  choice  of  a  controller  gain  to  satisfy 
the  Lyapunov  equations,  given  the  filter  gains.  The  sensitivity-shaped  filter  is  guar¬ 
anteed  to  be  asymptotically  stable  by  this  method,  as  long  as  the  controller  gains  are 
chosen  accordingly.  The  problem  is  posed  in  this  manner  to  evaluate  the  estimators 
relatively  independently  of  the  controllers.  Section  4.5  contains  a  closed-loop  system 
analysis  using  controller  gains  chosen  via  the  Lyapunov  method  of  Section  4.4.  In 
section  4.6,  a  sensitivity  analysis  with  respect  to  variations  in  the  structural  frequen¬ 
cies  up  to  20  %  is  included  for  each  type  of  estimator.  The  performance  criterir.  used 
is  the  root  mean  square  (rms)  of  the  error.  The  suboptimal  filters  are  shown  to  be 
less  sensitive  in  the  position  estimates.  Chapter  5  summarizes  the  comparison  of  the 
performance  of  various  estimator  types  under  perfect  and  imperfect  model  informa¬ 
tion  conditions.  Suggestions  for  future  research  in  the  area  of  estimator  types  applied 
to  LFSS  are  offered. 
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Chapter  2 


Opitmal  And  Suboptimal 
Estimation  Theory 


2.1  Introduction  to  Estimator  Types 


Since  not  all  desired  parameters  of  a  LFSS  are  known,  it  is  necessary  to  estimate  them. 
An  introduction  to  different  types  of  estimators  is  presented.  Imagine  a  generic  sys¬ 
tem  in  the  real  world,  call  it  system  G.  Suppose  it  is  composed  of  many  different 
components  with  generic  parameters,  call  them  P,  all  varying  in  time  and  space. 
To  be  able  to  predict  anything  about  system  G,  a  mathematic  model  Gm  represent¬ 
ing  the  physics  behind  it  must  be  formulated.  Since  the  actual  modelling  of  such  a 
system  may  be  difficult,  if  not  impossible,  many  simplifications  are  made  to  enable 
realistic  implementation  into  a  computer.  A  mathematical  set  of  equations  now  exists 
representing  system  G  and  its  motion  with  respect  to  some  reference  point  of  view 
as  it  moves  through  time  and  space.  Since  the  model  Gm  is  simplified,  eventually 
the  dynamics  of  Gm  will  depart  from  the  actual  dynamics  of  system  G;  this  may 
be  corrected  using  some  type  of  measuring  devices  and  feeding  the  measurements  to 
Gm.  Perhaps  parameters  P  are  desired  that  are  not  directly  measurable  quantities. 
Perhaps  more  accurate  information  is  needed  of  parameters  P  in  between  the  mea¬ 
surement  updates.  It  is  for  these  reasons  that  estimation  of  system  G  parameters  P 
is  introduced. 
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The  system  G  can  often  be  represented  as  a  finite-dimensional  linear  time-invariant 
system  in  state  space  form  as 


x(t)  =  Ax(t)  +  Bu(t)  (2.1) 

y(t)  =  Cx(t)  (2.2) 

with  x  the  state  vector,  u  the  input  vector,  and  A,  B  and  C  the  constant  plant,  control 
and  measurement  matrices,  respectively.  Here,  the  parameters  P  to  be  estimated  axe 
the  states  x  =  [ii,X2,  •••  xn]T •  To  estimate  the  states,  a  filter  is  formulated  as 


x(t)  =  AFx(t)  +  BFu(t)  +  I<F(y(t)  ~  (2.3) 

W)  =  Crm{t)  (2.4) 

where  *  is  the  estimated  state  vector,  AF,  Bp  and  CF  are  the  constant  estimator 
plant,  control  and  measurement  matrices,  respectively,  y  is  the  estimated  measure¬ 
ment  vector,  and  KF  is  the  filter  gain.  KF  may  be,  in  general,  time-varying  but  will 
be  assumed  time-invariant  for  this  study. 

The  problem  at  hand  is  the  choice  of  KF  which  will  determine  the  state  estimate  x 
and,  hence,  the  performance  of  the  filter.  It  is  proposed  to  investigate  several  methods 
of  determining  I(F.  The  optimal  filter,  in  the  sense  of  minimizing  the  mean  square 
value  of  the  error  e  =  x  —  x  for  a  given  set  of  known  model  parameters,  is  given  by 
Kalman  [6j.  This  is  a  full-order  centralized  filter  whereby  all  information  is  processed 
centrally  in  the  optimal  manner  stated  above.  Quite  often,  a  full-order  Kalman  filter 
is  too  cumbersome  and  costly  to  implement  on  a  computer.  A  reduced-order  Kalman 
filter  is  much  more  practical  in  this  case. 
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Suppose  now  that  central  processing  is  not  desirable  or  possible,  or  that  com¬ 
putation  costs  are  to  be  minimized.  Suppose  system  G  exists  in  such  a  way  as  to 
form  neat,  concise  subsets  of  the  states,  all  combining  to  form  the  total  system  G. 
A  separate  filter  could  be  designed  for  each  subset  of  system  G  which  would  process 
information  locally  -  information  could  then  be  contained  within  each  subset  with¬ 
out  the  need  to  feed  information  over  the  entire  system  G.  Such  a  set  of  filters  is 
a  decentralized  filter.  Advantages  of  a  decentralized  filter  include  enforcement  of  a 
desired  architecture,  and  usually  a  reduction  in  order  of  the  filters  needed  to  represent 
each  subsystem.  Because  of  this  reduction  in  order,  a  decentralized  filter  is  easier  to 
design  and  implement;  reduction  in  cost  is  an  important  factor  here.  Of  course,  the 
decentralized  filter  will  not  be  optimal  in  the  sense  of  minimizing  the  mean  square 
error  if  compared  to  the  centralized  filter,  but  the  above  advantages  may  justify  the 
degradation  of  performance. 

The  time-invariant  filters  described  above  are  categorized  as  optimal  (centralized 
full-order  filters)  or  suboptimal  (reduced-order  filters  and  decentralized  filters).  Also 
to  be  considered  are  filters  where  Kp  is  not  the  Kalman  gain,  but  rather  is  a  gain 
used  to  shape  filter  responses  or  some  other  set  of  specified  criteria. 


2.2  Centralized  Estimation  Theory 


The  full-order  centralized  estimation  equations  are  given  by  equations  (2.3)  and  (2.4) 
where  x  =  (zi,x2,...,  xn]T  and  n  is  the  order  of  the  modelled  system  Gm.  The 
optimal  gain,  in  the  sense  of  minimizing  the  mean  square  error,  given  a  set  of  model 
parameters,  is  the  Kalman  gain  [6] 


Kf  =  PCpR-1 


(2.5) 


where  P  is  the  solution  to  the  Riccati  equation 
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*  *  m  a  rrt  *  ,  a  a 

AFP  +  PAj.  -  PCJ-R-'CfP  +  q  =  0 


(2.6) 


and  Q  and  R  are  the  covariance  matrices  of  the  state  and  measurement  noises.  Equa¬ 
tion  (2.6)  assumes  a  steady  state  filter  exists,  which  implies  [ AF,BF ]  and  [Af,Cf]  are 
controllable  and  observable  pairs,  respectively.  An  example  of  a  centralized  system  is 
given  in  Figure  2.1. 

The  reduced-order  Kalman  filter  is  similar  to  the  full-order  case  with  correspond¬ 
ing  reduced-order  matrices 


xr  =  AFtxt  +  BFtut  +  I<F(yr  -  yr)  (2.7) 

yr  =  CFrx  r  (2.8) 

where  xr  is  the  reduced-order  state  estimate  vector  and  AFr,  BFt  and  CFr  are  the 
reduced-order  filter  plant,  control  and  measurement  matrices,  respectively.  The  input 
and  the  measurement  vectors  may  or  may  not  be  reduced-order,  depending  on  the 
existing  substructures  -  it  may  be  desirable  to  control  a  certain  part  of  a  LFSS  using 
only  local  measurements  and  actuators.  The  reduced-order  matrices  are  formed  by 
retaining  only  the  estimated  states  which  make  up  the  reduced-order  filter.  It  is  up 
to  the  discretion  of  the  designer  to  decide  which  states  are  most  important  to  retain. 
Many  techniques  of  model  reduction  applied  to  LFSS  exist  [15-20]  and  will  not  be 
discussed  in  this  study.  KF  is  the  reduced-order  Kalman  gain  determined  by 


I<F  =  PrCFrR;1  (2.9) 

where  Pr  is  the  solution  to  equation  (2.6)  with  the  appropriate  reduced-order  matrices 
substituted.  Qr  and  Rr  are  the  reduced-order  covariance  matrices.  This  reduced-order 
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filter  is  now  a  suboptimal  filter.  Error  has  been  introduced  because  certain  states  of 
the  system  have  been  neglected.  While  the  filter  is  no  longer  optimal,  it  is  still  the 
“best”  filter  for  a  system  of  the  same  order  as  the  reduced-order  model,  in  the  sense 
of  minimizing  the  mean  square  error. 

Suppose  Kf  is  now  chosen  to  shape  the  filter  response  of  equations  (2.7)  and  (2.8). 
It  may  be  advantageous  to  have  the  magnitude  of  the  filter  frequency  response  peak  in 
a  given  frequency  range.  It  may  be  desired  to  have  the  phase  of  the  frequency  response 
lie  in  a  certain  region.  This  Kp  then  has  been  chosen  to  “shape”  the  frequency 
response  of  the  filter  in  some  way  and  will  be  hereafter  named  a  frequency-shaped 
gain.  The  resulting  filter  is  a  suboptimal  filter  since  it  is  not  optimal  for  a  given  set 
of  model  parameters. 

2.3  Decentralized  Estimation  Theory 

It  is  assumed  that  the  system  dynamics  can  be  reduced  to  p  decoupled  subsets  of  nk 
equations  for  k  —  1,  ...,p,  where  no  further  reduction  within  each  of  the  p  subsets 
is  possible.  Note  that  the  sum  of  nk  over  k  is  n,  the  order  of  the  system.  Each 
decentralized  filter  is  formed  from  some  combination  of  the  p  subsets  where  the  equa¬ 
tions  for  any  two  of  the  decentralized  filters  axe  not  necessarily  decoupled.  The  ith 
decentralized  filter  is  given  by 

x,  =  A{Xi  +  BiUi  +  KfXVx  ~  Vi)  (2.10) 

and 

Vi  =  CiXx  (2.11) 

The  ith  estimated  state  vector  is  X;,  the  ith  plant,  control,  and  measurement  matrices 
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are  A;,  Bi,  and  C{  respectively,  and  the  ith  filter  gain  is  Kpt-  The  states  of  the 
decentralized  filter  subsets  may  be  chosen  using  model  reduction  techniques.  If  the 
decentralization  is  spatial,  different  modes  may  be  important  at  different  locations 
along  the  LFSS.  Or,  decentralization  may  be  used  to  reduce  the  overall  number  of  filter 
computations.  Each  filter  may  be  full-order  but,  by  taking  advantage  of  properties 
which  usually  exist  with  the  subset  architecture,  each  filter  will  most  probably  be 
reduced-order.  The  Kpi  may  be  chosen  to  be  reduced-order  Kalman  filter  gains  or 
some  criteria-shaped  gains.  These  filters  are  suboptimal  due  to  decentralization, 
reduction  in  order,  and,  if  criteria-shaped  gains  are  used,  not  minimizing  the  mean 
square  error.  An  example  of  a  decentralized  system  is  shown  in  Figure  2.2. 
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Figure  2.2:  Decentralized  filter 
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Chapter  3 


Modal  Analysis  And  Perfect 
Information  Results 


3.1  Modal  Analysis  Theory 

To  implement  a  control  system  for  a  flexible  structure,  a  dynamic  model  of  the  struc¬ 
ture  must  first  be  formed.  It  is  advantageous  to  express  the  model  in  modal  form, 
as  modal  control  -  the  damping  of  selected  modes  of  vibration  -  is  often  used.  The 
basic  equations  of  motion  for  a  structure  are 


[M]q  +  [Cl  j  +  \K\q  =  Q 


(3.1) 


where  [A/],  [C]  and  [K]  are  the  mass,  damping  and  stiffness  matrices,  respectively, 
q  is  the  displacement  vector  of  the  structure,  and  Q  is  the  generalized  force  vector 
acting  upon  the  structure. 

The  equations  of  motion  in  the  form  of  (3.1)  are  generally  coupled.  Implementing 
the  method  of  modal  analysis  as  given  in  [5]  involves  introducing  the  linear  transfor¬ 
mation 
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«(<)  =  ww) 


(3.2) 


where  [U]  is  the  modal  matrix  and  77(f)  is  a  generalized  coordinate.  Equation  (3.1) 
then  becomes 


[M)[U]f,  +  [C][U]i 7  +  [I<][U]V  =  Q.  (3.3) 


At  this  point,  a  justifiable  assumption  is  made  that  [C\  «  0.  A  term  will  be 
added  at  a  later  time  to  represent  the  small  inherent  damping  in  the  LFSS  without 
destroying  the  decoupled  form  of  the  equations.  Then 


wmv + mv\i  =  q 


(3.4) 


Premultiplying  equation  (3.4)  by  [U]T  yields 


[U]T{M][U)f,  +  [U]t[K)[U}ti  =  [U]tQ.  (3.5) 


The  amplitude  of  the  modal  matrix  [U]  is  arbitrary.  To  simplify  equation  (3.5), 
introduce  the  normalizations 


{V\t{M\[U)  =  [/] 


(3.6) 


[t/]T[/f][C/)  =  dtaglJf) 


(3.7) 


14 


where  [M]  and  [A']  are  positive-definite  symmetric  matrices,  and  u>,  axe  the  eigenval¬ 
ues.  Letting  [U^Q  =  N  results  in 


rj  +  diag(ui)r)  =  N 


(3.8) 


or 


Vi{t)  +“iV{t)  =  Ni(t) 


(3.9) 


for  the  ith  mode.  A  numerical  example  of  normalizing  [U]  will  be  given  later. 

As  a  representation  of  a  large  flexible  structure,  a  pinned-pinned  beam  is  consid¬ 
ered,  as  shown  in  Figure  3.1.  The  equation  for  an  Euler  beam  is 


at2 


dt* 


(3.10) 


where  m  is  the  mass  of  the  beam,  El  is  the  stiffness  of  the  beam,  y(x,  t)  is  the 
displacement  at  location  x  and  time  t  of  the  beam,  and  f(x,t)  is  the  force  distribution 
applied  to  the  beam.  To  solve  equation  (3.10),  an  assumption  that  free  vibration 
occurs,  i.e.  f(x,t)  =  0,  is  made.  If  y(x,t)  is  separable  in  time  and  space  as 


y{x,t)  =  q(t)U(x), 


(3.11) 


then  the  solutions  to  equation  (3.10)  obey 


q  -f-  u2q  =  0 


(3.12) 
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and 


t/(4)  +  {m^-)U  =  0.  (3.13) 

hi 

Solving  the  pinned-pinned  beam  problem  for  the  natural  frequencies,  w,,  and  the 
mode  shapes,  t/,-,  in  equations  (3.12)  and  (3.13)  yields  the  values 


u>i  =  (iT)2sqrt[ 


El 

(mL)‘ 


■],  i  —  1 » 2, ...,  oo 


(3.14) 


Ui(x)  =  Aisin(—). 


(3.15) 


The  mode  shapes  U,  are  arbitrary  in  amplitude  A{.  Normalizing  such  that 


(3.16) 


gives 

Ui{x)  =  sqrt[—]sin(~).  (3.17) 

The  total  displacement  and  velocity  of  the  beam  at  any  point  can  be  represented  as 
the  sum  of  the  mode  shapes  multiplied  by  the  generalized  coordinates: 


y{x,t)  =  J2Ui(x)Vi(t)  (3-18) 

t=i 
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(3.19) 


i=l 


The  dynamic  model  of  the  beam  has  now  been  expressed  in  modal  form  where  r},(t) 
and  r]i(t)  are  the  modal  positions  and  velocities  of  the  ith  mode,  respectively. 


3.2  Problem  Formulation 


The  simple  beam  example  used  throughout  this  study  is  the  pinned-  pinned  beam 
with  an  impulse  force  applied,  as  shown  in  Figure  3.1.  Velocity  sensors  and  collocated 
actuators  and  sensors  are  investigated  here  to  maintain  positive  realness  of  the  plant, 
as  discussed  in  Chapter  4. 

Let  the  state  *  be  defined  as 


*  =  [rjum  ■  ■  ... :  Vn,r)n}T 


(3.20) 


where  rji  and  t),  are  the  ith  modal  amplitude  and  velocity,  respectively.  The  state 
space  formulation  then  becomes 


x  =  Ax  4-  Bu 


(3.21) 


y  =  Cx 


(3.22) 


with 
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A  =  diag(Ai)  where  Ai  = 


0  1 

-w?  -20, 


(3.23) 


B  =  [0,  Ci 


o,cnf 


(3.24) 


and 


C  =  [0,Cx  :...:0,C2].  (3.25) 

Note  that  the  matrix  C  takes  the  above  form  for  a  velocity  sensor  and,  since  actuators 
and  sensors  are  collocated,  B  =  CT.  The  damping  term  £;  is  included  in  the  dynamics 
matrix  A  here  to  represent  the  slight  damping  inherent  in  the  structure  without 
destroying  the  decoupled  form  of  the  equations.  The  structural  model  is  now  in  state 
space  form  with  the  modal  positions  and  velocities  as  the  state  elements. 

Numerical  values  for  the  structural  parameters  are  given  in  Appendix  A. 

3.3  Estimation  Performance  Results  for  Perfect 
Model  Information 


Now  that  a  preferred  model  of  the  structure  exists,  estimation  of  the  state  is  de¬ 
sired.  Since  modal  control  is  to  be  implemented,  the  modal  positions  and  velocities 
are  required.  These  parameters  are  not  readily  available  for  measurement  and  must 
be  estimated.  The  filter  types  to  be  considered  are  full  and  reduced-order  central¬ 
ized  estimators,  full  and  reduced-order  decentralized  estimators,  and  a  reduced-order 
decentralized  estimator  with  frequency-shaped  gains.  In  this  chapter,  the  estima¬ 
tor  performances  are  compared  assuming  structural  frequency  model  information  is 
known  exactly. 
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The  basic  centralized  estimator  is  given  by  equations  (2.3)  and  (2.4)  and  is  pre¬ 
sented  again  in  equations  (3.26)  and  (3.27): 


x  =  Afx  +  Bfu  +  Kp(y  -  y). 


(3.26) 


y  =  CFx 


(3.27) 


where 


x  =  :  ...  :  T)n,rjn]T 


(3.28) 


Ap  =  diag(Api)  where  AFi  = 


0  1 
-w?  -2 C,-Wi 


(3.29) 


Bf  =  [0,  Ci  : ... :  0,  Cn]T 


(3.30) 


and 

r> 


Cp  =  [0,  Ci  : ... :  0,0]- 


(3.31) 


The  optimal  time-invariant  gain,  in  the  sense  of  minimizing  the  mean  square  error, 
given  a  set  of  model  parameters,  is  the  Kalman  gain 


Kf  =  PCpR-1 


(3.32) 


where  P  is  the  solution  to  the  Riccati  equation 


AfP  +  PAI-PCJ'R-1CfP  +  Q  =  0  (3.33) 


and  Q  and  R  are  the  covariance  matrices  of  the  state  and  measurement  noises. 

Two  cases  with  identical  filters  but  different  system  models  are  considered.  The 
first  uses  noise  in  the  model  when  calculating  the  Kalman  gain,  even  though  no 
noise  is  present  in  the  system.  The  noise  used  when  obtaining  the  Kalman  gain 
is  called  pseudo-noise.  The  second  case  determines  the  Kalman  gain  assuming  noise 
covariances  that  correspond  to  actual  noise  in  the  system.  The  same  noise  covariances 
are  used  in  both  cases  and  numerical  values  are  presented  in  Appendix  A.  The  first 
filter  is  actually  not  an  optimal  filter  -  the  optimal  steady  state  Kalman  gains  for  a 
full-order  filter  with  no  noise  in  the  system  are  zero.  The  second  filter  is  the  optimal 
time-invariant  filter  for  the  system  with  the  modelled  noise.  The  statistical  noise 
information  for  each  mode  is  given  in  Appendix  A. 

The  results  of  the  centralized  estimation  are  shown,  in  conjunction  with  the  actual 
motion,  in  Figures  3.2  to  3.3.  The  first  six  modes  of  the  system  model  describe  the 
actual  motion  of  the  system  adequately.  However,  the  motion  is  dominated  by  the 
first  two  modes.  The  six-mode  centralized  estimator  is  full-order,  while  the  one-mode, 
two-mode  and  four-mode  estimators  are  reduced-order.  All  the  centralized  estimators 
use  the  Kalman  gains  (full  or  reduced-order)  as  the  filter  gains  and  are  termed  Linear 
Quadratic  (LQ)  estimators.  Note  that  the  centralized  estimate  follows  the  actual 
motion  closely  for  the  two  and  higher-mode  filters.  The  one-mode  filter  contains 
second  mode  and  higher  frequencies  in  the  error,  due  to  low  damping  and  to  the 
high  frequency  content  of  the  measurement  input  which  contains  all  of  the  modes. 
Since  the  beam  motion  is  dominated  by  the  first  two  modes,  this  phenomenon  does 
not  appear  in  the  two-mode  and  higher-mode  filters.  Examining  Figures  3.2  and 
3.3,  two  facts  of  value  are  noted.  The  velocity  estimates  appear  to  be  “better”  than 
the  position  estimates  -  this  is  because  velocity  sensors  are  used.  Also,  the  trends 
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appearing  in  the  system  without  noise  do  not  change  when  system  noise  is  included. 

For  the  decentralized  case,  the  full-order  decentralized  estimator  will  be  identical 
to  the  full-order  centralized  estimator  if  information  is  not  combined  in  a  suboptimal 
manner  and  if  the  Kalman  gain  is  chosen.  Hence,  the  results  for  this  filter  will  not  be 
presented.  As  an  arbitrary  reduced-order  decentralized  case,  the  identical  example  of 
Figure  3.1  is  used,  but  with  two  velocity  sensors  located  as  shown  in  Figure  3.4.  One 
sensor  provides  output  for  one  estimator  while  a  second  sensor  provides  output  for 
a  second  estimator.  The  first  estimator  models  and  estimates  only  the  first  mode  of 
the  system  while  the  second  estimator  models  and  estimates  only  the  second  mode 
of  the  system.  This  case  represents  the  worst  scenario  -  here,  the  dominant  modes 
of  the  motion  are  assumed  to  be  unknown.  Hence,  the  result  is  two  single-mode 
decentralized  estimators.  The  ith  estimator  is  given  as 


Xi  =  AFiii  +  BFtu{  +  KFi{Vi  -  Vi)  (3.34) 

y,  =  CFiXi  (3.35) 

with  BFi  =  [0 ,Ci]T,CFi  =  [0,  C,],  and  i,  =  [7,-,  r/,]T.  The  results  of  the  two  de¬ 
centralized  estimators  can  be  combined  through  summation  if  there  are  no  common 
states,  or  if  there  are  common  states,  through  any  desired  estimation  algorithm  [21] 
such  as  Mean  Square  Error  estimate  (MSE),  Maximum  Apriori  estimate  (MAP),  or 
Maximum  Likelihood  estimate  (MLE). 

Using  Kalman  gains  calculated  with  pseudo-noise  in  the  model  but  with  no  noise 
present  in  the  system,  Figure  3.5  shows  the  oscillatory  motion  of  the  decentralized 
estimates  about  the  actual  motion.  Examining  the  separated  first  and  second  mode 
estimates  of  the  decentralized  system,  shown  in  Figure  3.6,  it  can  clearly  be  seen 
that  the  error  between  the  decentralized  estimate  and  the  actual  motion  oscillates  at 
undesirably  high  frequencies.  This  is  due  to  low  damping  and  to  the  high  frequency 
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measurement  input  -  the  same  phenomenon  that  occured  with  the  one-mode  cen¬ 
tralized  filter.  Since  this  oscillation  causes  large  errors  in  the  position  and  velocity 
estimates,  it  is  advantageous  to  remove  the  high  frequency  error  content  from  the 
decentralized  estimate  through  the  shaping  of  the  Bode  magnitude  plots  of  the  error 
transfer  functions,  called  error  amplitude  shaping. 

The  error  transfer  functions  for  the  decentralized  estimates  are  derived  from  the 
relation 

e  =  Hx  -  x  (3.36) 


or 


E(s)  =  [si  -  AFcL\-l[HA  -  AFclH  -  I<FC]X{s)  (3.37) 

where  H  is  a  matrix  of  ones  and  zeros  that  selects  the  elements  of  the  reduced-order 
state  from  the  actual  state,  and  AFcl  =  AF  —  KFCF  is  the  closed- loop  filter  dynamics 
matrix.  By  varying  the  gains  in  equation  (3.37),  the  frequency  response  of  the  filters 
can  be  shaped  to  the  desired  response. 

For  the  one-mode  decentralized  filter,  the  denominator  of  the  transfer  matrix  is  a 
second-order  polynomial  given  by 

s2  +  (2£i oji  -f  KiC\)s  +  w2(l  —  K\C\).  (3.38) 


Comparing  equation  (3.38)  with  the  general  equation  for  a  second-order  system 


s 2  +  (2£u;n)s  +u>l 


(3.39) 
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yields  the  natural  frequency  and  damping  ratio  of  the  filter  response  as 

uin  =  u>\  sqrt(l  —  KiCi)  (3.40) 

and 

_  2(1q;1  +  A"2Ci 

^  “  2u>n 

The  designer  can  therefore  choose  the  desired  natural  frequency  and  damping  ratio 
of  the  filter  via  the  gains  I\\  and  A'2.  The  transfer  functions  to  be  considered  are 
the  error  in  the  first  mode  estimated  position  and  velocity  due  to  the  actual  second 
mode  velocity,  and  the  error  in  the  second  mode  estimated  position  and  velocity  due 
to  the  actual  first  mode  velocity.  The  Bode  magnitude  plots  of  these  error  transfer 
functions,  shown  in  Figure  3.7,  indicate  second  mode  and  higher  frequencies  in  the 
first  and  second  mode  filters.  By  varying  the  gains  in  (3.18),  the  frequency  response  of 
the  error  transfer  functions  may  be  shaped  as  shown  in  Figure  3.8.  Figure  3.9  shows 
the  separated  first  and  second  estimated  modes  of  motion  with  frequency-shaped 
filters.  It  can  be  seen  that  much  of  the  high  frequency  error  has  now  been  removed. 
The  total  actual  motion  of  the  system  and  decentralized  estimation  results  with  both 
filters  frequency-shaped  are  shown  in  Figure  3.10.  Much  improvement  over  Figure  3.5 
is  observed. 

Note  that  the  previous  decentralized  results  correspond  to  a  system  with  no  noise 
present.  Now  that  a  set  of  satisfactory  frequency-shaped  gains  has  been  chosen  for 
the  decentralized  filter,  these  same  gains  are  used  in  a  system  with  noise  as  given  in 
Figure  3.3.  The  results  are  shown  in  Figure  3.11. 


(3.41) 
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Position  at  .3L  (ft)  Position  at  .3L  (ft) 


Figure  3.2a:  Centralized  position  estimation  results  for  the  system  with  no  noise 
(pseudo  noise  used  to  determine  gains)  and  perfect  model  information 


25 


Velocity  at  .3L  (ft/sec)  Velocity  at  .3L  (ft/sec 
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Figure  3.3a:  Centralized  position  estimation  results  for  system  with  noise  and  perfect 
information 
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Velocity 


Figure  3.4:  Pinned-pinned  beam  with  impulse  F  applied  ;  Sensor  1  -  Connected  with 
estimator  for  1st  mode  ;  Sensor  2  -  Connected  with  estimator  for  2nd  mode 
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Figure  3.6:  Decentralized  estimates  using  LQ  gains  of  the  1st  and  2nd  modes  of 
motion  for  system  with  no  noise  and  perfect  information 
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Frequency  w  (rad/sec) 


Figure  3.8a:  Bode  magnitude  plots  for  error  in  first  mode  estimate  due  to  second  mode 
velocity.  System  has  no  noise  and  perfect  model  information.  Filter  is  decentralized 
filter  with  shaped  gains. 


Velocity  at  .3L  (ft/sec)  Position  at  .3L  (ft) 


Figure  3.9.  Shaped  decentralized  estimates  of  the  1st  and  2nd  modes  of  motion  for 
system  with  no  noise  and  perfect  information 
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Velocity  at  .3L  (ft/sec)  Position  at  .3L  (ft) 


2 


Figure  3.10:  Total  shaped  decentralized  estimates  for  system  with  no  noise  and  perfect 
information 
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Chapter  4 

Robustness  And  Model  Error 
Sensitivity  Analysis 


4.1  Sensitivity-Shaped  Gains  Selection 


The  previous  analysis  has  been  done  for  full  and  reduced-order  linear  quadratic  (LQ) 
estimators  and  for  a  reduced-order  decentralized  estimator.  No  precautions  have  been 
taken  to  a^ure  these  filters  give  good  performance  in  the  presence  of  modelling  errors. 
A  sensitivity-shaped  filter  is  now  formulated  to  address  this  problem.  A  set  of  filter 
gains  must  be  determined  to  ensure  the  system  guaranteed  asymptotic  stability  in 
the  presence  of  modelling  errors. 

The  modelling  errors  are  chosen  as  errors  in  the  structural  frequencies  u>  of  up 
to  20  %■  This  is  felt  to  be  representative  of  the  actual  uncertainties  of  the  first  few 
structural  modes  of  a  LFSS. 

There  is  probably  an  inexhaustible  supply  of  methods  which  will  determine  a  sat¬ 
isfactory  set  of  gains;  this  study  will  employ  one  method.  The  method  is  based  upon 
minimizing  the  sensitivity  of  the  filter  to  variations  in  the  structural  frequencies.  The 
purpose  of  this  method  is  to  improve  filter  performance  when  imperfect  knowledge  of 
the  system  frequencies  exists.  While  structural  mode  shape  information  is  also  uncer- 
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tain,  it  is  assumed  in  this  study  that  mode  shapes  are  known  exactly.  The  equations 
of  motion,  from  equations  (2.1)  to  (2.4),  in  state  variable  form  are 


x  =  Ax  +  Bu  (4.1) 

y  =  Cx  (4.2) 

where  A,  B  and  C  are  the  system  plant,  control  and  measurement  matrices.  The 
filter  is  represented  by 


x  =  Apx  +  Bfu  +  I<f(v  ~  y)  (4-3) 

y  =  Cpx  (4.4) 

with  Ap,  Bp  and  Cp  as  the  filter  plant,  control  and  measurement  matrices  and  I\p 
as  the  filter  gain  matrix.  The  actual  system  plant  matrix  A  will  contain  the  actual 
frequencies  u>i,  which  may  vary  from  the  modelled  frequencies  by  20  %.  The  filter 
plant  matrix  Ap  will  contain  the  modelled  frequencies  u><.  It  is  interesting,  then, 
to  note  the  transfer  matrix  between  the  actual  state  x  and  the  estimated  state  *. 
Since  this  relationship  is  Multi-Input  Multi-Output  (MIMO),  the  transfer  matrix  is 
composed  of  the  individual  transfer  functions  between  the  inputs  and  outputs.  If 
the  estimator  was  of  the  same  order  as  the  system  plant,  if  noise  was  not  taken  into 
account,  and  if  there  were  no  modelling  errors,  then  this  transfer  matrix  should  be 
the  identity  matrix.  It  is  proposed  to  place  emphasis  on  this  point  and  to  attempt  to 
design  a  filter  with  this  property. 

Let  Hp(s)  represent  the  filter  transfer  matrix  between  the  filter  output,  x,  and  the 
filter  input,  the  state  measurement  y.  Assuming  that  the  number  of  actuators  and 


40 


sensors  are  identical,  or  dim(u )  =  dim(y),  an  expression  for  HF(s)  can  be  determined 
by  defining  the  control  vector  u  as 


u  =  [C(sl  -  A)-1i?]-1y 


(4.5) 


and  substituting,  with  equation  (4.4),  into  (4.3)  to  obtain  equation  (4.6).  The  result 
is 


Hf(s)  =  [si  -Af  +  I<fCf]~1[Bf{C(sI  -  A)_1B}_1  +  I<F).  (4.6) 


The  system  remains  open-loop,  but  the  feedback  loop  is  driven  by  the  same  reference 
input  as  the  system  plant.  A  block  diagram  of  this  open-loop  system  is  shown  in 
Figure  4.1.  Recall  that  the  plant  matrices  A  and  AF  are  partitioned  into  matrices 
with  Ai  and  AFi,  respectively,  along  the  diagonals  and  zeros  on  the  off-diagonals 
resulting  in 


A  =  diag{Ai)  where  A\  = 


0  1 
-(S?  -2(iW, 


(4.7) 


and 


Af  =  diag(AF' )  where  AFi  = 


0  1 
-u>7  -2C.U4 


(4.8) 


Through  matrix  partitioning,  the  plant  and  filter  equations  can  be  written  as 


Xi  =  A{X  +  BiUi 


(4.9) 
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Vi  =  CiXi 


(4.10) 


x=AFtXi  +  BFlu{  +  I<F,(yt  ~  Vi)  (4-H) 

Vi  =  CFiX{  (4.12) 

with  Xi  =  [ T}i,i};]T,B ,  =  Bf>  =  [0 ,Ci]T,Ci  =  CFt  =  [0,  C,j,  and  x,  =  [t),,?),]7’.  The 
system  control  and  measurement  matrices  B  and  C  are  affected  by  variations  in  the 
mode  shapes,  but  not  the  structural  frequencies.  Hence,  Ui  terms  do  not  appear  in 
these  matrices.  In  this  study,  it  is  assumed  that  mode  shapes  are  perfectly  known. 

The  transfer  matrix  S(s )  between  the  estimated  state  x  and  the  actual  state  x  is 
just 


S(s)  =  Hf(s)C. 


(4.13) 


It  is  desired  to  force  this  matrix  S(s)  as  close  to  the  identity  matrix  as  possible. 
First,  let  it  be  noted  that  5(s)  will  only  be  square  if  the  system  model  has  the  same 
number  of  modes  representing  it  as  the  filter  model.  Since  the  extra  terms  rendering 
the  matrix  non-square  could  just  be  forced  to  zero,  the  same  number  of  modes  will 
be  used  for  simplicity.  Starting  with  a  one-mode  representation  of  the  system  plant 
and  filter,  S(s)  has  the  elements 


5(s) 


m/vi  ji/vi 
iih\  vJm  . 


(4.14) 


These  elements  can  be  analytically  determined  from  equations  (4.6)  and  (4.13)  as 
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mhi 


■s^  +  (2CiWi  +  K^C^s  -f-  w2(l  —  K\C) 
~  52  +  (2Ciwt  +  I<2C)3  +  w*(l  -  K\C) 


(4.15) 


mlv  i  = 


_ l(i(KiC  —  l)(^i  —&i) _ 

s2  +  (2CiWi  +  K2C)s  +  w2(l  -  I<\C) 


(4.16) 


VilVi  = 


_ s(^i  ~  u\) _ 

S2  +  (2ClWl  +  K2C)s  +  U>2(1  -  KXC) 


(4.17) 


*  ..  s2  +  (2Ciu>i  +  K3C)s  +  a;2(l  -  KXC) 
VllVl  s2  +  (2CiWj  +  I<2C)a  +  w2(l  -  IUC) 


(4.18) 


where  Kp ,  =  [K\,  K-Jf  and  C\  =  (0,  C]T.  Note  that  if  d>i  =  u>\,  then  equations  (4.15) 
to  (4.18)  reduce  to 


mhi  =  1 


(4.19) 


771/771  =  0 


(4.20) 


Vilvi  =0 


(4.21) 


vjii  =  1 


(4.22) 


In  an  attempt  to  drive  the  off-diagonal  terms  of  equations  (4.15)  to  (4.18)  to  zero  and 

the  diagonal  terms  to  one,  the  value  of  K\  =  7;  is  chosen.  Equations  (4.15)  to  (4.18) 

C 

become 
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fi\h\  =  i 


(4.23) 


7)1/771  =  0 


(4.24) 


‘  /  -  *(^1  ~  ^1) 

Vl,Tl1  s(s  +  2CiWi  +  I<2C) 


(4.25) 


5(5  -4-  2^i  wi  +  X2C) 

s(5  +  2^iU7i  +  I<2C) 

Choosing  a  value  for  I<2  which  shifts  the  pole  (— 2Ci<*>i  +  K2C)  sufficiently  to  the  left 
will  ensure  that  the  response  due  to  that  pole  dies  out  quickly.  Equations  (4.23)  to 
(4.26)  then  approach  those  of  (4.19)  to  (4.22). 

So  far,  the  method  used  above  seems  to  be  rather  specific.  It  assumes  the  number 
of  modelled  system  plant  modes  is  equal  to  the  number  of  modelled  filter  modes. 
It  assumes  that  only  one  mode,  the  same  mode  for  the  system  plant  and  filter,  is 
modelled.  In  the  2-mode  case,  using  the  first  two  modes  as  an  example,  the  values  of 
the  filter  gains  are  much  more  difficult  to  ascertain.  No  easy  choices  are  immediately 
seen  from  the  equations.  With  the  aid  of  a  symbolic  manipulator,  the  transfer  matrix 
between  x  and  x  is 


■  J 

I< 

L 

M  ' 

N 

0 

P 

Q 

R 

S 

T 

u 

V 

w 

X 

Y 

(4.27) 


where 
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J  —  S*  +  [2ClWi  +  2^2U>2  4"  C2K4  +  Ci/C2]33  + 

[4ClU>lC2w2  +  2Cl  +  2/C4C2C1W!  +  W|  {1  —  Ci/^i}  +  w|{l  —  C^/G}]^2  4- 

[Ci-^2wf  4"  C2K4OJ1  +  2£2u>2u)2{1  —  C\K])  4-  2£iu/ju>|{1  —  C2/^3}]s  4- 
tt>2u>2{l  —  C\K\  —  C2i^3}.  (4.28) 


K  —  [u>2  —  wjjs3  4-  [{u;2  —  u2}{C2/C4  4-2C2w2}]s2  4-  [^{l  ~  C2^3}{w2  —  ^i)]5  (4.29) 

L  —  [Cii^3{u;2  —  <Z>2}]s2  4-  [{w2  —  u>i }{2Ci/^3^2<^2  4*  (4.30) 

M  =  [CXI<<{ w}  -  *2}]s2  +  [Ct^3u;2{w2  -  a;2}]*  (4.31) 


^  —  [2Ci{l  ~  CiKi}{ui  —  wi}]32  4-  [{u>i  —  Wi}{4CiC2a,2(l  —  CXI<X)  4-  2C2/^4Ci}]s  4- 
2{u>i  -  wtHCiwftl  -  CXI<X  -  C2K3)}  (4.32) 


0  —  s4  4"  [2^i ujx  4"  2^2^2  4"  C2K4  4-  Cj/^l^3  4* 

[4CiwiC2w2  4-  2CxI\2^2'jJ7  4-  2K4C2C1&1  4-  w2{l  —  CXKX}  4-  w2{l  —  C2/i3}]s2  4- 

\CxK2^j\  4-  C2K4U1 2  4-  2£2w2u;J{1  —  CXKX}  4"  2(i<Dic<;2{1  —  C2/^3}].s  4- 

u>xu>2  { 1  —  CXKX  —  C2K3}.  (4.33) 


P  —  [2Ci/v3^i{a;i  —  wx}]s2  4-  [{wj  —  i^\}2^xCx{2K3^2uj3  4-  /^JJs  (4-34) 
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Q  —  [2Ci/y4£i{u;i  —  wi}]s2  —  [2CiK3^i4{ui  —  u>i}]s 


(4.35) 


R  —  [C\Ki{oJ2  —  4"  [{^2  —  4  (4.36) 

s  =  [C2I<2{ul  -u\}}s2  +  [C2K\<jj\{u\  -  <4 }}s  (4.37) 


•s 4  +  +  2£2ui2  -+•  C2K±  4  C\K2]s3  + 

[4Cia,iC2a;2  4-  2ClK2(^2u>2  4  2K4C2(\Ui  +  o?2{1  —  C\K\ }  -4  —  C2K$}\s2  + 

\CiK2u\  -4  C2K\u\  4-  2^2^2wi{1  —  C\K\ }  -4  2Ci\uj\Hi\  { 1  —  C2Kz}\s  4- 
44  —  C2Ki  —  C2K3)'  (4.38) 


4  ~4W  +  {{4  ~4}{CiI<2  4  4iu>i}]s2  4  [u;2{l  -  C^H*}  -  4}}s  (4.39) 


V  —  \1C2K\(12{u2  —  o>2}]52  4  [{^2  —  w2}2(12C2{2K\(i\uji  +  A2}]5  (4.40) 


W  —  [2C2K2C2&2  —  W2}]s2  4  \{y>2  ~  (4-41) 


[2C2 { 1  —  C2IU){u2  —  a>2}]s2  4  [{^2  —  W2}{4CiC2Wi(1  —  C2K3)  4  2CiK2C2}]s  4 
2{u2  -  W2}{C2W?(1  -  C,I<2  -  C2/f3)}  (4.42) 


Y  =  s4  4-  -f  2(2^2  +  C2K4  4-  CiK^s3  4- 

[4^x011^2^2  4*  20x1^2^2^2  4"  2K4C2 Ci^i  +  ^1  { 1  —  C\K\}  -f  o)j {1  —  C2/C3}]52  + 

\C\K2<jj\  4-  C2K4U1  +  2^2^2wi{1  —  C\K\ }  +  2(>iijj\uj\ { 1  —  C2/V3U3  4- 

u\u\ { 1  -  Cx  Kx  -  C2IU }  •  (4.43) 


=  54  +  [2£iU>i  +  2^2^2  +  C2K4  4*  C\K^S3  + 

[4^xWx^2w2  4*  2Cx R 2(2^2  4-  2K4C2C1W1  4-o?j{1  —  C\K\ }  4-  <*>2 {1  ”  C2K3W31  4- 


[Ci/fjWj  4"  C2K4U1  4*  2^2W2Wx{1  —  CxKx}  4-  2£xu>iu;2{1  —  O2K2WS  4- 
u\u\{l  -  CxKx  -  C2K3}- 


(4.44) 


Trying  to  force  S  to  the  identity  matrix  is  not  as  obvious  as  in  the  one-mode 
case.  However,  choosing  K\  =  —  and  K%  —  30  as  in  the  one-mode  case,  and 

1  1 

K$  —  Ka  =  70,  yields  a  two- mode  sensitivity-shaped  filter  which  produces  more 
C2 

accurate  velocity  estimates  than  the  one-mode  sensitivity-shaped  filter.  Since  the 
pinned-pinned  beam  is  dominated  by  the  first  two  modes  of  motion,  it  is  not  necessary 
to  produce  a  filter  of  three  modes  or  higher.  However,  it  is  conjectured  that  the  gains 
for  the  n-mode  sensitivity-shaped  filter  can  be  determined  via 


*»-■  -  cl 


K21  >  w,+l 


(4.45) 


for  i  =  1,  ...,n. 


Using  the  six-mode  model  for  the  system  plant,  as  outlined  in  Chapter  3,  the 
time  response  of  the  sensitivity-shaped  estimators,  the  central  LQ  estimators,  and 


the  decentral  frequency-shaped  estimator  are  compared.  The  estimators  all  contain 
the  modelled  frequencies  a?,.  Systems  with  and  without  noise  present  are  included. 
For  the  case  with  no  system  noise,  given  in  Figure  4.2,  the  sensitivity-shaped  filters 
estimate  position  more  accurately,  in  the  root  mean  square  (rms)  error  sense,  than 
the  LQ  filters  with  pseudo-noise  or  the  decentral  filter.  The  fact  that  the  sensitivity¬ 
shaped  filters  perform  better  than  the  full-order  LQ  filter  is  due  to  the  fact  that 
pseudo-noise  was  introduced  to  obtain  the  LQ  gains.  The  optimal  time-invariant 
filter  of  zero  gains  is  just  the  model  with  nominal  frequencies,  and  is  also  included  in 
Figure  4.2.  The  velocity  estimates  of  the  one-mode  and  two-mode  sensitivity-shaped 
filters  are  comparable  to  the  two-mode  and  four-mode  LQ  filters,  respectively.  For 
the  system  including  noise,  given  in  Figure  4.3,  the  one-mode  sensitivity-shaped  filter 
performs  comparable  to  the  four-mode  and  two-mode  LQ  filters  for  position  and 
velocity  estimates,  respectively.  Increasing  the  sensitivity-shaped  gain  k2  above  the 
given  value  will  yield  even  better  velocity  estimates.  Thus,  the  one-mode  sensitivity¬ 
shaped  filter  is  at  least  preferred  over  a  filter  of  two  modes  -  a  saving  in  computations 
is  achieved  here.  The  two-mode  sensitivity-shaped  filter  estimates  position  more 
accurately  than  any  other  filter  investigated.  The  fact  that  the  sensitivity-shaped 
position  rms  error  values  are  lower  than  those  of  the  optimal  full-order  filter  is  due  to 
the  fact  that  the  optimal  filter  is  time-invariant  and  also  to  the  statistical  anomaly  in 
generating  noise  in  the  simulation.  The  two-mode  sensitivity-shaped  filter  accuracy 
of  the  velocity  estimates  lies  between  the  two-mode  and  four-mode  LQ  filters.  It  is 
important  to  note  that  these  results  are  independent  of  sensor  position.  Regardless 
of  the  sensor  position,  the  same  trends  in  performance  of  the  filters  are  recognized. 

4.2  Positive  Real  Systems 

A  set  of  filter  gains  to  decrease  filter  frequency  sensitivity  has  been  chosen  on  the  basis 
of  estimation  performance.  Is  the  resulting  filter  going  to  yield  a  robust  closed-loop 
system?  Before  analyzing  the  robustness  of  the  system,  it  is  important  to  introduce 
positivity  concepts.  The  following  Definitions  1  and  2  and  Theorems  1  and  2  are 
presented  in  [4]. 
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Definition  1:  A  square  transfer  matrix  H(s)  is  called  positive  real  if 

(a)  H(s)  has  real  elements  for  real  s,  and 

(b)  H(s)  has  elements  which  are  analytic  for  3fJ[s]  >  0,  and 

(c)  H*(s)  +  H(s )  is  positive  semidefinite  for  9£[s]  >  0,  where  *  denotes  the 
complex  conjugate  transpose. 


Definition  2:  A  square  transfer  matrix  H(s)  is  strictly  positive  real  if 

(a)  H(s)  has  real  elements  for  real  s,  and 

(b)  H(s)  has  elements  that  are  analytic  for  3?[s]  >  0,  and 

(c)  H*(ju>)  +  H(ju})  is  positive  definite  for  all  real  u>. 


Theorem  1: 

Given  the  square  transfer  matrices  G(s)  and  H(s)  in  the  feedback  system 
of  Figure  4.4,  the  system  is  asymptotically  stable  if  at  least  one  of  the 
transfer  matrices  is  positive  real  and  the  other  transfer  matrix  is  strictly 
positive  real.  If  both  transfer  matrices  G(s)  and  H (s)  are  only  positive 
real,  then  the  feedback  system  will  be  at  least  marginally  stable. 


In  order  to  satisfy  the  sufficient  condition  of  Theorem  1,  it  must  be  determined 
whether  G(s)  or  H(s)  is  strictly  positive  real  -  the  other  will  be  constrained  to  be 
positive  real.  Normally,  conditions  ensuring  a  strictly  positive  real  transfer  matrix 
G(s)  are  too  restrictive  on  the  model.  Hence,  G(s)  will  be  considered  the  positive 
real  transfer  matrix.  G(s)  is  not,  however,  positive  real  in  general  and  the  system 
must  be  modelled  in  a  certain  form  to  attain  this. 


Theorem  2: 

For  the  system  described  in  state  variable  form  in  equations  (4.1)  and 
(4.2),  the  transfer  matrix  G(s)  =  C(sl  -  A)~lB  is  positive  real  if  and 
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only  if  there  exists  a  symmetric  positive  definite  matrix  P  and  a  symmetric 
positive  semi-definite  matrix  Q  such  that 


pa  +  atp  =  -g 

(4.46) 

CT  =  PB. 

(4.47) 

The  equations  of  motion  of  the  LSS  in  terms  of  modal  variables  is 


T)  +  [2Cu>]j)  -1-  [w2]t7  =  N. 


(4.48) 


Equation  (4.48)  can  be  written  for  n  modes  in  the  state  space  form,  with 


*  =  [Vi,Vi  •  ••• :  Vn,Vn\ 


(4.49) 


and 


B  =  [0,  Ci  :  ... :  0,  Cn]T .  (4.50) 

Now,  from  Theorem  2,  for  a  positive  real  transfer  matrix  G(s),  conditions  (4.46)  and 
(4.47)  must  hold.  Because  of  the  form  of  the  A  matrix  in  equation  (4.7),  it  can  be 
shown  that  (4.46)  will  be  satisfied  for 


P  =  diag(Pi)  with  Pi  = 


(4.51) 


Q  =  diag(Qi)  with  Qi  = 


0  0 

0  4C,u;t 


(4.52) 
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Substitution  of  (4.50)  and  (4.51)  into  (4.47)  yields 


ct  =  b. 


(4.53) 


The  matrix  B  defined  in  (4.50)  operates  on  the  control  vector  u  to  indicate  how 
the  input  signals  act  on  the  system  states.  It  holds  the  information  on  the  actuator 
locations  and  orientations  through  the  generalized  force  vector  N.  The  matrix  C 
operates  on  the  state  vector  to  give  the  outputs.  It  will  therefore  hold  the  information 
on  the  sensor  locations  and  orientations.  Therefore,  using  CT  =  B  implies  that  the 
system  has  collocated  actuators  and  sensors,  and  the  above  analysis  indicates  that  this 
system  configuration  guarantees  a  positive  real  system  transfer  matrix  G(s).  Note 
that  the  form  of  C  given  by  (4.50)  and  (4.53)  restricts  the  system  sensors  to  rate 
sensors  only.  A  positive  real  G(s )  may  exist  for  non-collocated  actuators  and  sensors 
and/or  sensors  other  than  rate  sensors.  However,  these  limitations  have  been  adhered 
to  here  since  they  guarantee  G(s)  positive  real. 


4.3  Phase-shaping  Method  to  Test  Robustness 


Theorem  1  is  a  sufficient  condition  for  guaranteed  asymptotic  stability  in  the  presence 
of  modelling  errors.  For  collocated  actuators  and  rate  sensors,  the  system  transfer 
matrix  G(s)  is  positive  real.  Attempts  will  now  focus  on  finding  a  feedback  transfer 
matrix  H(s)  which  is  strictly  positive  real  as  given  by  Definition  2.  This  can  undoubt¬ 
edly  be  examined  in  numerous  ways.  One  approach  is  termed  the  “phase-shaping” 
method,  which  is  outlined  below. 

A  sufficient  condition  for  H(s)  to  be  strictly  positive  real  is  given  in  Theorem  3. 


Theorem  3: 

If  the  phase  of  each  transfer  function  element  of  the  filter  transfer  matrix 
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Hf(s)  lies  between  ±90  degrees,  and  if  the  controller  gains  I\  are  positive, 
then  H(s)  =  KHf(s)  is  strictly  positive  real.  Note  that  if  the  endpoints 
±90  degrees  are  included,  then  H(s)  is  positive  real.  A  proof  of  this 
theorem  is  contained  within  Appendix  B. 


The  focus  of  this  study  is  on  the  performance  of  the  optimal  and  suboptimal 
filters.  The  open-loop  case  will  be  presented  here,  as  well  as  the  closed-loop  case  with 
generic  controllers.  It  is  assumed  that  the  state-feedback  controller  gain  K  must  be 
chosen  to  ensure  that  Theorem  3  is  satisfied.  In  this  way,  filter  performance  can  be 
evaluated  as  independently  as  possible  of  the  controller  used. 

Referring  to  Figure  4.5,  the  feedback  transfer  matrix  H(s)  is  the  combination  of 
the  estimator  transfer  matrix,  Hp(s),  and  the  controller  gain  matrix,  K.  H(s)  is 
written  as 


H{s)  =  KHf{s).  (4.54) 

The  dimensions  of  these  matrices  are  as  follows:  K  is  the  number  of  measurements 
by  twice  the  number  of  filter  plant  modelled  modes;  HF{s)  is  twice  the  number  of 
filter  plant  modelled  modes  by  the  number  of  measurements;  and  H(s)  is  the  number 
of  measurements  by  the  number  of  measurements. 

In  the  case  of  the  one-mode  sensitivity-shaped  filter  of  Section  4.1,  there  is  a 
phase  shift  of  -90  degrees.  This  indicates  positive  realness  rather  than  strict  positive 
realness.  The  Bode  phase  plots  for  the  transfer  function  elements  of  the  transfer 
matrix  between  the  filter  output  (estimated  state)  and  the  filter  input  (measurement) 
are  given  in  Figures  4.6  to  4.7. 

While  these  results  do  not  guarantee  asymptotic  stability,  two  things  must  be 
remembered.  For  the  one-mode  sensitivity-shaped  fiber,  Ff(.s)  is  guaranteed  to  be 
positive  real  by  the  sufficient  condition  of  the  phase-shaping  method  -  this  in  turn 
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guarantees  the  closed-loop  system  at  least  marginal  stability  by  the  sufficient  condi¬ 
tion  of  Theorem  1.  In  other  words,  a  conservative  constraint  has  been  imposed  on  top 
of  another  conservative  constraint.  The  phase-shaping  method  guarantees  that  the 
system  using  the  sensitivity-shaped  filter  is  at  least  marginally  stable.  These  results 
rely  upon  the  choice  of  K  to  satisfy  Theorem  3. 


4.4  Lyapunov  Method  to  Test  Robustness 


The  phase-shaping  method,  based  on  the  sufficient  condition  of  Theorem  3,  does  not 
indicate  whether  or  not  the  system  is  guaranteed  asymptotic  stability  in  the  presence 
of  modelling  errors.  Another  approach  is  to  investigate  Theorem  2  as  a  necessary  and 
sufficient  condition.  A  form  of  Theorem  2  is  given  by 


PAd  +  A%P  =  -(Q  -  I(TC  -  CTI()  (4.55) 


Ad  =  (A  -  BK) 


(4.56) 


I\T  =  PKf 


(4.57) 


where  P  is  a  symmetric  positive  definite  matrix  and  Q  is  a  symmetric  positive  semidef- 
inite  matrix.  For  the  sensitivity-shaped  gains  of  Section  4.1,  H(s )  is  strictly  positive 
real  if  a  K,P  and  Q  can  be  found  to  satisfy  equations  (4.55)  to  (4.57). 

Looking  at  the  one-mode  sensitivity-shaped  filter,  the  matrices  become 


P  = 


Pll  P 12 
Pi  2  P22 


(4.58) 
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q  _  9n  9i2 
6  9l2  ?22 


(4.59) 


Kf  =  [kfu  kf2]T  (4.60) 

I<  =  [kuk2],  (4.61) 

The  constraints  are 

9n  =  ~2pi2(u>i  +  fciCi)  (4.62) 

9i2  ==  k\Ci  +  Pu  —  Pi2(2CiWj  +  k2C\)  —  p22{(jj\  +  fciCi)  (4.63) 

922  =  2  fc2Cx  +  2(pi2  —  P22[2Ci^i  +  k2CV])  (4.64) 

k\  —  Pukfi  +  p22kf2  (4.65) 

k\  =  pi2kf\  +  p22kf2  (4.66) 

P  >  0  (4.67) 

Q>  0.  (4.6S) 


1  .4893 

.4893  13.99 


yields  the  controller 


Choosing  P  -  |  _  0038  .0522  J  and  ^  ~ 

gain  matrix  K  =  [.1108, 1.5654]  which  will  satisfy  equations  (4.55)  to  (4.57).  For  the 
set  of  sensitivity-shaped  gains,  then,  it  is  possible  to  choose  P,  Q,  and  I\  to  ensure 
H(s)  is  strictly  positive  real.  In  actuality,  there  are  many  different  choices  of  P,  Q, 
and  K  which  will  ensure  this,  and  Section  4.5  deals  with  this  topic. 


4.5  Closed-Loop  System  Analysis 


A  brief  analysis  of  the  effect  of  P  and  Q  on  the  controller  gain  matrix  K  and,  hence, 
the  closed-loop  system  performance,  is  presented.  Six  cases  axe  chosen  to  represent 
possible  variations  in  kx  and  fc2  and  are  listed  below  as  Cases  1  through  6  ir.  equa¬ 
tions  (4.69)  to  (4.74). 


1 

-.0038  ' 

,Q  = 

l 

.49 

-.0038 

.05 

.49 

13.95 

[  .1108  .5 

1 

1  .49 
.49  13.95 

I<  =  [  .1108  1.5654  ] 

1  .49 
.49  14.73 

J 

I<  =  [  .1108  24  ] 


’  1  .16  ' 
.16  .051 

,Q  = 

1 

.65 

.65 

13.98 

[  5  1.5654 

1  .49 

.49  14.1 


P  = 


1  -.0038 

.0038  .166 


,<?  = 


P  = 


1  -.0038 

.0038  .8 


P  = 


1 

.0038 


-.0038 

.05 


,Q  = 


(4.69) 


(4.70) 


(4.71) 


(4.72) 
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K  =  [55]  (4.73) 

1  .545 

545  13.99 

K  =  [  3.37  1.5654  ]  (4.74) 


P  = 


1  .105 

.105  .05 


The  closed-loop  time  response  for  each  of  these  cases  is  given  for  the  system  with  no 
noise  present  in  Figures  4.S  to  4.9.  Note  that  the  larger  the  control  forces,  the  faster 
the  actual  positions  and  velocities  are  damped,  as  is  expected.  Referring  back  to 
Section  4.3,  do  the  controller  gains  chosen  via  the  Lyapunov  method  also  satisfy  the 
controller  gain  requirements  of  the  phase-shaping  method?  Since  all  of  the  controller 
gains  are  positive,  then  Theorem  3  guarantees  H(s )  positive  real  which,  in  turn, 
guarantees  marginal  stability  by  the  phase-shaping  method  of  Section  4.3. 


4.6  Sensitivity  to  Model  Error  Analysis 


The  system  using  the  one-mode  sensitivity-shaped  filter  of  Section  4.1  and  the  con¬ 
troller  gains  of  Section  4.4  has  been  shown  to  be  robust  for  structural  frequency 
variations  of  20  %.  How  does  the  performance  of  this  filter  compare  with  those  using 
LQ  estimation  theory?  This  is  related  to  the  sensitivity  of  the  estimates  to  changes 
in  the  model  parameters,  namely,  the  structural  frequencies.  It  is  stated  in  [7]  that 
optimal  filters  are  more  sensitive  to  model  errors,  as  shown  in  Figure  4.10. 

The  sensitivity-shaped  filters  of  Section  4.1,  the  LQ  one,  two,  four  and  six-mode 
filters,  and  the  reduced-order  decentralized  filter  with  frequency-shaped  gains  are 
considered  here.  The  structural  frequencies  vary  from  the  modelled  frequencies  by 
±  20  %  in  5  %  increments.  Figure  4.11  shows  the  sensitivity  to  variations  in  the 
frequencies  u>  of  each  of  these  filters  in  the  position  and  velocity  estimates  with  no 
noise  modelled  in  the  system.  Model  reduction  among  the  LQ  filters  indicates  that  the 
lower  the  order  of  the  filter,  the  higher  the  rms  errors,  an>  -he  less  sensitive  to  changes 
in  u>.  This  corresponds  Lo  the  trends  shown  in  Figure  4.10.  For  the  position  estimates. 
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the  sensitivity-shaped  filters  perform  the  best  under  structural  frequency  variations. 
Once  again,  the  one-mode  and  two-mode  sensitivity-shaped  results  are  comparable 
to  the  two- mode  and  four-mode  LQ  filters,  respectively,  for  the  velocity  estimates. 
Figure  4.12  shows  the  sensitivity  to  variations  in  u>  of  the  sensitivity-shaped,  LQ  and 
decentralized  filters  for  the  system  with  modelled  noise  included.  The  sensitivity¬ 
shaped  filters’  position  estimates  are  less  sensitive  to  frequency  variations.  Although 
the  velocity  estimate  sensitivity  plot  no  longer  contains  a  minimum,  due  to  statistical 
anomalies  in  the  noise,  the  same  relative  performace  of  the  filters  is  noted,  and  the 
one-mode  and  two-mode  sensitivity-shaped  filters  again  perform  comparable  to  the 
two-mode  and  four-mode  LQ  filters,  respectively. 

In  general,  the  accuracy  of  the  position  estimates  of  the  sensitivity-shaped  filters 
are  superior  to  the  velocity  estimates.  This  is  because  the  entries  in  the  S  matrix 
of  equations  (4.14)  and  (4.27)  corresponding  to  position  are  identically  one  and  zero 
while  the  entries  in  the  S  matrix  corresponding  to  velocity  are  determined  by  how 
much  the  gain  k2i  shifts  the  corresponding  pole  to  the  left,  causing  its  effects  to  die 
out  at  different  rates.  The  gain  &2«  should  be  chosen  large  enough  to  yield  acceptable 
velocity  estimates;  if  there  is  a  bound  on  fc2l,  then  the  velocity  estimate  performance 
is  defined  via  this  bound. 
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Position  at  .3L  (ft)  Position  at  .3L  (ft) 


Figure  4.2a:  Position  estimation  results  for  various  filters  with  no  noise  in  the  system 
and  structural  frequency  variations  of  20  % 
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Position  at  .3L  (ft)  Position  at  .3L  (ft) 


Figure  4.2b:  Position  estimation  results  for  various  filters  with  no  noise  in  the  system 
and  structural  frequency  variations  of  20  % 
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Velocity  at  .3L  'ft/sec)  Velocity  at  .3L  (ft/sec 


30 


r 


Figure  4.2c:  Velocity  estimation  results  for  various  filters  with  no  noise  in  the  system 
and  structural  frequency  variations  of  20  % 


61 


Actual  6-mode  motion 
1-mode  LQ  estimate 
Decentral  freq. -shaped  2-mode 
estimate 


Actual  6-mode  motion 

1- mode  sens -shaped  estimate 

2- mode  sens -shaped  estimate 


Time  (sec) 

Figure  4. 2d:  Velocity  estimation  results  for  various  filters  with  no  noise  in  the  system 
and  structural  frequency  variations  of  20  % 


Figure  4.3a:  Position  estimation  results  for  various  filters  with  modelled  noise  included 
and  structural  frequency  variations  of  20  % 


Velocity  at  .3L  (ft/sec)  Velocity  at  .3L  (ft/sec) 


Figure  4.3b:  Velocity  estimation  results  for  various  filters  with  modelled  noise  in¬ 
cluded  and  structural  frequency  variations  of  20  % 
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Figure  4.3c:  Estimation  results  for  the  1- 
included  with  structural  frequency  variation 
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Figure  4.4:  General  feedback  control  concept 
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4.5:  General  feedback  control  concept  with  estimator  and  controller  in  feedback 
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Figure  4.6:  Bode  phase  diagram  of  outp 
mate)  to  inpu1  '  neasurement  y).  System 
variations  of  20  % 


Position 


Case  1 
Case  2 
Case  3 

(increasing  Case  #  indicates 
increasing  k  for  fixed  k, ) 


Time  (sec) 

Figure  4.8a:  Closed-loop  system  position  response  using  various  controllers  and  the 
1-mode  sensitivity-shaped  estimator  with  no  noise  in  the  system  and  structural  fre¬ 
quency  variations  of  20  % 


Velocity  at  .3L  (ft/sec) 


Time  (sec) 

Figure  4.Sb:  Closed-loop  system  velocity  response  using  various  controllers  and  the 
l -mode  sensitivity-shaped  estimator  with  no  noise  in  the  system  and  structural  fre¬ 
quency  variations  of  20  % 
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Position  at  .3L  (ft) 


Time  (sec) 

Figure  4.9a:  Closed-loop  system  position  response  using  various  controllers  and  the 
1-mode  sensitivity-shaped  estimator  with  no  noise  in  the  system  and  structural  fre¬ 
quency  variations  of  20  % 


Velocity 


Time  (sec) 

Figure  4.9b:  Closed-loop  system  velocity  response  using  various  controllers  and  the 
1-mode  sensitivity-shaped  estimator  with  no  noise  in  the  system  and  structural  fre¬ 
quency  variations  of  20  % 
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Figure  4.9c:  Closed-loop  system  position  response  using  various  controllers  and  the 
1-mode  sensitivity-shaped  estimator  with  no  noise  in  the  system  and  structured  fre¬ 
quency  variations  of  20  % 


Velocity  at  .3L  (ft/sec 


Figure  4.9d:  Closed-loop  system  velocity  response  using  various  controllers  and  the 
1-mode  sensitivity-shaped  estimator  with  no  noise  in  the  system  and  structural  fre¬ 
quency  variations  of  20  % 
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EXPECTED 


—  RANGE  of  — 

uncertasnmty 


Fi0ure  4.10.  Conceptual  example  of  designing  for  minimum  sensitivity  (crtrue  held 
constant,  crde3tgn  is  varied  y  K  tTUt 
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Position  rms 


Figure  4.11a:  Sensitivity  of  various  estimators  to  variations  in  structural  frequencies 
up  to  20  %.  System  does  not  contain  noise. 
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Velocity  rms 


Actual/ modelled  frequency 

°f  Vari0USKlimatMS  ‘°  variations  in  structural  frequencies 
up  to  20  Yo.  System  does  not  contain  noise. 
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Figure  4.12b:  Sensitivity  of  various  estimators  to  variations  in  structural  frequencies 
up  to  20  %.  Noise  present  in  system. 
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Chapter  5 


Conclusions  and 
Recommendations 


Full  and  reduced-order  centralized  and  decentralized  estimators  applied  to  a  flexible 
structure  have  been  evaluated.  When  perfect  model  information  exists,  the  full  order 
Kalman  filter  is  the  optimal  filter  in  the  sense  of  minimizing  the  mean  square  value  of 
the  error.  While  errors  are  introduced  via  model  order  reduction  and  decentralization, 
these  filters  generally  require  less  computation  and  design  costs.  A  more  extensive 
investigation  into  the  errors  introduced  through  decentralization,  including  a  wider 
range  of  decentralized  estimators,  is  suggested. 

In  actuality,  perfect  model  information  does  not  exist  and  the  Kalman  filter  is  no 
longer  optimal  under  these  modelling  uncertainties.  A  sensitivity-shaped  filter  has 
been  formulated  to  reduce  sensitivities  to  variations  in  the  structural  frequencies  and 
its  estimates  are  compared  to  the  centralized  Linear  Quadratic  (LQ)  filters  and  the 
decentralized  frequency-shaped  filter,  all  with  20  %  structural  frequency  errors.  The 
sensitivity-shaped  filters’  position  estimates  are  superior  to  all  other  filters  examined 
when  modelling  errors  exist.  The  one-mode  and  two-mode  sensitivity-shaped  filters’ 
velocity  estimates  are  comparable  to  the  two-mode  and  four-mode  centralized  LQ 
filters,  respectively.  Future  work  should  include  an  investigation  of  the  proposed 
general  method  to  calculate  gains  of  an  n-mode  filter. 
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Robustness  of  the  closed-loop  system  with  the  one-mode  sensitivity-shaped  filter 
has  been  evaluated  using  a  phase-shaping  method  which  guaranteed  the  closed-loop 
system  marginal  stability,  and  using  the  Lyapunov  method  which  guaranteed  the 
closed-loop  system  asymptotic  stability,  both  under  structural  frequency  errors  of  20 
%.  Both  methods  also  assumed  the  controller  gains  met  certain  given  requirements. 
Several  controllers  were  investigated  to  examine  the  closed-loop  system  response  - 
their  gains  were  chosen  via  the  Lyapunov  method.  The  controller  gains  chosen 
via  the  Lyapunov  method  also  satisfied  the  phase-shaping  method  requirements  for 
controller  gains.  Preliminary  studies  indicate  a  need  for  further  research  into  possible 
controller  types  and  resulting  closed-loop  system  responses. 

A  sensitivity  to  structural  frequencies  analysis  has  been  conducted  for  all  filter 
types  introduced.  Suboptimal  filters  are  shown  to  be  less  sensitive  to  structural  fre¬ 
quency  variations  than  optimal  filters.  The  sensitivity-shaped  filters’  position  rms 
errors  were  lower  than  the  centralized  LQ  and  decentralized  frequency-shaped  filters 
when  modelling  errors  existed.  The  one-mode  and  two-mode  sensitivity-shaped  fil¬ 
ters’  velocity  rms  errors  were  comparable  to  the  two-mode  and  four-mode  centralized 
filters,  respectively,  which  indicates  a  saving  in  computation  through  order  reduction. 

Finally,  the  estimators  evaluated  in  this  study  have  been  applied  to  a  simple 
pinned-pinned  beam.  It  is  important  to  extend  these  results  to  a  structure  which 
would  more  closely  represent  a  LFSS,  such  as  a  truss. 


« 
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Appendix  A 


mode  1 

mode  2 

mode  3 

mode  4 

mode  5 

mode  6 

frequency 

7.0689 

28.2757 

63.6202 

113.1026 

176.7228 

254.4809 

damping  ratio 

.01 

.02 

.03 

.04 

.05 

.06 

measurement 

[0,  4.4501: 

0,  5.2314: 

0,  1.6998: 

0,  -3.2332: 

0,  -5.5007: 

0,  -3.2332] 

matrix 
pseudo  noise 


Q  =  .5  *  I(order  2 n),  R  =  .2  *  I(order  dim(i/))  ,  n  =  #  of  modes 


mean  of  state  noise  (0.0) 

-.011  .009  -.008  .015  .014  -.035  -.032  .026  -.004  -.081  -.014  .009 


variance  of  state  noise  (.5) 

.499  .510  .530  .502  .500  .431  .483  .501  .470  .449  .465  .491 


t 


standard  deviation  of  state  noise  (.707) 

.706  .714  .728  .708  .707  .657  .695  .708  .685  .670  .682  .701 
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mean  of  measurement  noise  (0.0) 
-.019 


variance  of  measurement  noise  (.2) 
.217 


* 


standard  deviation  of  measurement  noise  (.447) 
.466 


Filter  types 

Filter  gains 

1-mode  LQ 

[-1.372,  5.741]t 

2-mode  LQ 

[-.460,  10.521,  -1.502,  12.616]T 

4-mode  LQ 

[-.212,  10.829,  -.784,  37.037,  -.869,  61.831,  1.490,  -8.485]: 

6-mode  LQ 

[-.116,  10.778,  -.406,  40.362,  -.267,  70.596,  .446, 
-127.575,  1.235,  -109.889,  1.118,  132.774]T 

Decentral  frequency-shaped  2-mode 
mode  1:  [-.01,  .01]T 

mode  2:  [-.01,  .1] 

T 

1-mode  sensitivity-shaped 
[.2247, 30] r 


2-mode  sensitivity-shaped 
[.2247, 30,. 1912,  70]t 


» 


< 
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Appendix  B 


* 


t 


Proof  of  Theorem  3: 

Hf(s)  is,  in  the  case  of  the  one-mode  sensitivity-shaped  filter,  a  column  matrix 
with  transfer  function  elements  [Hp1  Hf3].  If  the  phase  of  a  transfer  function  element 
lies  between  ±  90  degrees  not  including  the  endpoints,  then  the  transfer  function  is 
a  strictly  positive  real  function.  Also  for  this  case,  the  controller  gain  matrix  is  given 
by  [ki  k2\-  If  and  k2  are  both  positive,  then  the  feedback  transfer  matrix  H(s)  = 
KHf(s)  is  the  sum  of  two  strictly  positive  real  transfer  functions,  kxHp2  +  k2Hp2. 
The  statement  to  prove,  then,  is  that  the  sum  of  strictly  positive  real  functions  is  a 
strictly  positive  real  function. 


Let  H(s)  =  Hi(s)  +  H2(s)  +  ...  +  Hr{s).  Then 

(a)  If  Hi(s),H2(s)1  ...,Hr(s)  have  real  elements  for  real  s,  then  H(s)  also 
has  real  elements  for  real  s. 

(b)  If  Hx(s),  ...,Hr(s)  have  elements  analytic  for  $R(s)  >  0,  then  H(s)  has 
real  elements  for  3£(s)  >  0. 

(c)  By  the  definition  of  a  positive  real  function,  for  all  x  ^  0,  and  for  all 
real  u> , 

xm[Hi(ju)  +  H;(ju)]x  >  0,  ...,x*[Hr(jw)  +  H;(ju)]x  >  0. 

Then  the  sum  of  these  functions  is  also  positive 
x‘[Hi{ju)  +  H\{ju)\x  +  ...  -(-  x*[Hr{ju)  +  H*{ju)]x  >  0. 

This  can  be  factored  as 
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x  [Hi(ju)  +  ...  +  Hr(ju)  +  Hi(ju)  +  ...  +  H*(jui)]x  >  0  or 
x'[II(ju)+ H'(ju)\x  >  0  which  states  that  H(s )  is  a  positive  real  function. 


» 


T 


* 

< 


*U.S.Government  Printing  Office:  1990— 748-056/24470 


88 


